The Las Campanas Infrared Survey. III. The i7-band Imaging 
Survey and the Near-Infrared and Optical Photometric Catalogs 
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ABSTRACT 

The Las Campanas Infrared Survey, based on broad-band optical and near- 
infrared photometry, is designed to robustly identify a statistically significant and 
representative sample of evolved galaxies at redshifts z > 1. We have completed 
an H-hand imaging survey over 1.1 square degrees of sky in six separate fields. 
The average 5 a detection limit in a four arcsecond diameter aperture is ~ 
20.8. Here we describe the design of the survey, the observation strategies, data 
reduction techniques, and object identification procedures. We present sample 
near-infrared and optical photometric catalogs for objects identified in two survey 
fields. The optical images of the Hubble Deep Field South region obtained from 
the literature reach 5 a detection thresholds in a four arcsecond diameter aperture 
of [/ ~ 24.6, B ~ 26.1, V ~ 25.6, R ~ 25.1, and / ~ 24.2 magnitude. The 
optical images of the Chandra Deep Field South region obtained from our own 
observations reach 5 a detection thresholds in a four arcsecond diameter aperture 
of V ~ 26.8, R ~ 26.2, / ~ 25.3, and z' ~ 23.7 mag. Wc perform object detection 
in all bandpasses and identify ^ 54,000 galaxies over 1,408 arcmin^ of sky in the 
two fields. Of these galaxies, ~ 14,000 are detected in the iJ-band and ~ 2,000 
have the colors of evolved galaxies, I — H ^ 3, z ^ 1. We find that (1) the 
differential number counts N{m) for the if-band detected objects has a slope 
of dlog N{m)/dm = 0.45 ± 0.01 mag-^ at // ^ 19 and 0.27 ± 0.01 mag-^ at 
H 1^ with a mean surface density ~ 7, 200 degree"^ mag~^ ai H = 19. In 
addition, wc find that (2) the differential number counts for the if-band detected 
red objects has a very steep slope, d\og N{m;I — H ^ 3) /dm = 0.84 ± 0.06 
mag~^ at if ^ 20 and 0.32 ± 0.07 mag~^ a.t H ^ 20, with a mean surface density 
~ 3, 000 degree"^ mag~^ at H — 20. Finally, we find that (3) galaxies with red 
optical to near-IR colors {I — H > 3) constitute fa 20% of the if-band detected 
galaxies at H ^ 21, but only Ki 2% at H ^ 19. We show that red galaxies 
are strongly clustered, which results in a strong field to field variation in their 
surface density. Comparisons of observations and predictions based on various 
formation scenarios indicate that these red galaxies are consistent with mildly 
evolving early-type galaxies at 2; ~ 1, although with a significant amount of 
on-going star formation as indicated by the large scatter in their V — I colors. 

Subject headings: catalogs — cosmology: observations — galaxies: evolution — surveys 
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1. INTRODUCTION 

Different galaxy formation scenarios provide distinct predictions for tlie space densities 
and masses of evolved galaxies at redshifts beyond z = 1. In monolithic collapse scenarios, 
massive galaxies form early over a dynamical time (e.g. Eggen, Lynden-Bell, & Sandage 1962) 
and passively evolve to the present epoch. The co-moving space density of evolved galaxies 
is therefore expected to remain constant and the intrinsic luminosities of these galaxies are 
expected to increase gradually with increasing redshift. In hierarchical formation scenarios, 
massive galaxies form later through mergers of low-mass galaxies over a significant fraction 
of the Hubble time (White & Rees 1978). The bulk of star formation and mass assembly 
occurs much later and the co-moving space density of evolved galaxies is expected to decline 
steeply with redshift a.t z > 1. While various surveys have identified some luminous galaxy 
populations at z > 2 (e.g. Steidel et al. 1999; Blain et al. 1999), the connection between 
these galaxies and typical present-day galaxies is unclear. On the other hand, comparisons 
of statistical properties of evolved high-redshift galaxies and local eUiptical galaxies may 
provide a direct means of discriminating between competing galaxy formation scenarios (e.g. 
Kauffmann & Chariot 1998). 

Evolved galaxies may be characterized by their intrinsically red colors due to the lack 
of ongoing star formation that provides most of the UV light in typical galaxy spectral 
energy distributions. At redshifts beyond unity, where the UV region of the spectrum is 
redshifted to the visible, optical galaxy surveys are insensitive to evolved galaxies, because 
of their steep UV spectra. Identifications of evolved high-redshift galaxies therefore rely 
on observations carried out at near-infrared wavelengths. Various studies based on near- 
infrared surveys have, however, yielded inconsistent measurements of the space density of 
evolved galaxies (Kauffmann, Chariot, & White 1996; Totani & Yoshii 1997; Zepf 1997; 
Franceschini et al. 1998; Benitez et al. 1999; Menanteau et al. 1999; Barger et al. 1999; 
Schade et al. 1999; Broadhurst & Bouwens 2000; Daddi et al. 2000; Martini 2001a). Because 
most of the existing deep, near-infrared surveys observe a sky area of only a few tens of 
square arc-minutes, the discrepancy is likely due to significant field-to-field variation (e.g. 
Daddi et al. 2000; McCarthy et al. 2001a,b). Furthermore, because evolved galaxies lack 
prominent narrow-band spectral features in the UV spectral range, this discrepancy is also 
likely due to selection biases (e.g. Totani & Yoshii 1997). Additional uncertainty may also 
arise due to the presence of dusty star forming galaxies that exhibit red colors resembling 
evolved high-redshift galaxies (e.g. Graham & Dey 1996; Small, Ivison, & Blain 1997). 

A complete survey of evolved galaxies over a large sky area is needed to address these 
issues. Deep, wide-field near-infrared surveys have only recently become feasible due to 
the advent of large format near-infrared cameras (Beckett et al. 1998; Elston 1998; Persson 
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et al. 2001). In addition, various groups have demonstrated that distant galaxies may be 
accurately and reliably identified using photometric redshift techniques that incorporate op- 
tical and near-infrared broad-band photometry (Connolly et al. 1997; Sawicki, Lin, & Yee 
1997; Lanzetta, Fernandez- Soto, & Yahil 1998; Fernandez-Soto, Lanzetta, & Yahil 1999; 
Fernandez-Soto et al. 2001). For galaxies lacking strong emission or absorption features, we 
may still be able to estimate redshifts based on the presence/absence of spectral discontinu- 
ities using photometric redshift techniques. Therefore, a wide-field optical and near-infrared 
imaging survey, with sufficiently accurate colors to yield reliable photometric redshifts, of- 
fers an opportunity to identify statistically significant and representative samples of evolved 
galaxies at 2; > 1. 

We initiated the Las Campanas Infrared (LCIR) Survey to obtain and analyze deep 
near-infrared images and complementary optical images in the V, R, I, and z' bandpasses 
over one square degree of sky at high galactic latitudes (Marzke et al. 1999; McCarthy et 
al. 2001a,b). The survey is designed to identify a large number of red galaxies at 2; > 1, 
while securing a uniform sample of galaxies of all types to 2; ~ 2 using broad-band optical 
and near-infrared colors. In particular, this program utilizes one of the largest near-infrared 
cameras available (ClRSl), which produces an image of 13' x 13' contiguous field of view in 
a sequence of four pointings. The LCIR survey may serve to bridge the gap between very 
wide field surveys (e.g. NOAO Deep- Wide; Jannuzi et al. 1998) and very deep, smaller field 
near-infrared imaging surveys (e.g. ESO Imaging Survey, da Costa 2000; NTT Deep Field, 
Moorwood, Cuby, & Lidman 1998). 

The primary objectives of the program are (1) to examine the nature of the red galaxy 
population and identify evolved galaxies at redshifts 2; > 1, (2) to study the space density 
and luminosity evolution of evolved galaxies at redshifts z ^ 2, and (3) to measure the 
spatial clustering of evolved high-redshift galaxies (McCarthy et al. 2001a,b; Firth et al. 
2001), thereby inferring merging rates of these galaxies for constraining theoretical models. 
The galaxy catalog together with an accompanied photometric redshift catalog will also be 
used (4) to assess luminosity and luminosity density evolution for galaxies of different types 
over the redshift range between 2; = to 2 (Chen et al. in preparation) and (5) to study the 
bright-end luminosity functions for galaxies and QSOs at redshifts beyond z = 4.5. Finally, 
galaxies identified in the survey will be compared with objects identified in a companion 
VLA survey (6) to study the spatial distribution of weak radio sources 10 /xJy; Yan et al. 
in preparation). 

To date we have completed an if-band imaging survey over 1.1 square degrees of sky 
to a mean 5 a detection limit in a four arcsecond diameter aperture of if ~ 20.8. Here we 
present sample optical and near-infrared photometric catalogs for galaxies identified in two 
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of our fields: the Hubble Deep Field South (HDFS) and Chandra Deep Field South (CDFS). 

These are the two fields for which a complete set of near-infrared and optical images arc 
available. Wc have identified ~ 24,000 galaxies over 847 arcmin^ in the HDFS region, of 
which ~ 6, 720 are detected in the if-band survey, and ~ 30, 000 galaxies over 561 arcmin^ 
in the CDFS region, of which ~ 7, 400 are detected in the /J-band survey. Among the 
if-detcctcd galaxies in the two regions, ~ 2, 000 have colors that match evolved galaxies, 
/ — ;^ 3, at redshifts z 1. A complete near- infrared and optical photometric catalog 
of the two fields may be accessed at http:/ /www.ociw.edu/lcirs/catalogs.html. Catalogs for 
the other survey fields will be released as they become available. 

We describe the survey design and observation strategies in § 2 and data reduction 
and image processing in § 3. The procedures for identifying objects in multi-bandpasses, 
estimates of the completeness hmits of the catalogs, and an assessment of the rehability of 
the photometric measurements are presented in § 4. The format of the optical and near- 
infrared photometric catalogs is described in § 5. Finally, we discuss statistical properties of 
the if-band selected galaxies in § 6. An independent analysis based on objects identified in 
the HDFS region is also presented in Firth ct al. (2001). Wc adopt the following cosmology: 
VIm = 0.3 and I^a = 0.7 with a dimensionless Hubble constant h = Ho/{lQQ km s~^ Mpc~^) 
throughout the paper. 



2. Observations 

2.1. Definition of the Las Campanas Infrared Survey 

The goals of the LCIR survey are to robustly identify a statistically significant and 
representative sample of evolved galaxies at redshifts z > 1, as well as to secure a uniform 
sample of galaxies of all types at intermediate redshifts based on broad-band optical and 
near-infrared colors. 

To determine the image depths necessary to reach our goals, we estimated the expected 
near-infrared brightness and optical and near-infrared colors of evolved high-redshift galaxies 
using a series of evolutionary models. We considered a non-evolving model using an empirical 
early-type galaxy template presented in Coleman, Wu, & Weedman (CWW, 1980) and two 
evolving models for galaxies of solar metallicity formed in a single burst of 1 Gyr duration at 
redshifts Zf = 5 and Zf = 10, respectively (Bruzual & Chariot 1993). We scaled these models 
to have Lh = LhJ at ^ = (see also Marzke et al. 1999). Figure 1 shows the predicted 



tWe adopt Mh, = 



-22.9 + Slog /i from De Propris et al. (1998). 
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redshift evolution of the observed i?-band magnitudes and optical and near-infrared colors 
of evolved galaxies based on the three models. Adopting Hq — 70km s"-*^ Mpc~^, wc found 
that in the absence of dust an evolved galaxy may have an apparent if-band magnitude 
between H = 19.5 and 20.3 and / — if ~ 3 at redshift 2; ~ 1, and an apparent if-band 
magnitude between H = 20.7 and 22.5 and an I — H color between ~ 4.5 and 5 at redshift 
2; ~ 2. The color threshold I — H ^ 3 therefore defines our selection criterion for evolved 
galaxies at redshifts z > 1, which according to Figure 1 is similar to selections made based 
on R — Kg ^ 5. To identify a large number of the evolved galaxy population, we therefore 
choose a target 5 a survey depth of if ~ 21, more than one magnitude fainter than an L* 
elliptical galaxy at 2; ~ 1. We also require that the optical imaging survey reach consistent 
depths based on the predicted optical and near-infrared colors. The actual depth that we 
achieved varied from field-to-field and often fell somewhat short of our design goals. 

To determine the survey area, we first estimated the number of evolved galaxies that are 
required to obtain a significant signal in the clustering analysis. At a projected co-moving 
correlation length Vp = 1 Mpc, wc found that a sample of ~ 3000 galaxies are needed to reach 
the 5(J level of significance (Marzke et al. 1999). The surface density of galaxies satisfying 
i — if 3 is, however, very uncertain. Different measurements based on K-hand surveys 
have been reported by several groups, ranging from 0.3 arcmin"^ to 3 arcmin"^ at Ks < 20 
(Cowie et al. 1994; Elston et al. 1997; Barger et al. 1998). Adopting a median surface density 
1 arcmin"^ and a mean near-infrared color if — fCg ~ 1 for evolved galaxies at 2; > 1, we 
concluded that one square degree of sky at a 5 cr limiting magnitude if ~ 21 will yield ~ 3000 
red high-redshift galaxies. 

The LCIR survey may serve to bridge the gap between the on-going NOAO Wide-Field 
Survey that images 18 square degrees of sky at a shallower depth (Jannuzi et al. 1998) and 
the existing deep, smaller field near-infrared imaging surveys such as the NTT Deep Field 
(Moorwood, Cuby, & Lidman 1998). We have selected six distinct equatorial and southern 
fields at high galactic latitudes distributed in right ascension. In Table 1, we list the fields, 
their J2000 coordinates, and galactic latitudes. 

2.2. Infrared Imaging Observations 

The if -band imaging survey was carried out at the du Pont 2.5 m telescope at Las Cam- 
panas using the Cambridge Infrared Survey Instrument (CIRSI). This instrument contains 

a sparse mosaic of four 1024 x 1024 Rockwell Hawaii HgCdTe arrays (Beckett et al. 1998: 
Persson et al. 2001). At the Casscgrain focus of the du Pont telescope, the plate scale of 
the camera is 0.199 arcsec pixel"^. Some of the observations were made with a doublet field 
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corrector in front of the dewar window, which produced a plate scale of 0.196 arcsec pixel"^. 
The spacing between each array is 90% of an array dimension. Therefore, one pass of four 
pointings ~ 192 arcsec apart in a square pattern is needed to fill in the spaces between the 
arrays. A filled mosaic covers ~ 13' x 13'. This defines our unit field area, henceforth referred 
to as a "tile". 

Observations of individual pointings were composed of nine sets of three to five exposures 
35 to 45 s in duration with dither offsets of between 8 and 13 arcsec in a rectangular pattern. 
For a typical H-hand sky brightness of 14.7 magnitudes, we estimated that a total exposure 
time of 80 minutes was needed to reach a 5 cr limiting magnitude of i/ ~ 20.3 in a four arcsec 
diameter aperture. The target integration time was achieved by making three passes on each 
field. Because it was often impossible to complete three passes of a field within one night, 
the observations of any given field area were spread over several nights. Standard stars from 
Persson et al. (1998) were observed three to five times per night. Flat field observations were 
made using a screen hung inside the dome and were formed from differences of equal length 
exposures with the dome lamps on and off. 

Most of our fields are composed of four contiguous tiles arranged to cover a roughly 
25' X 25' field. The HDFS field contains seven tiles arranged in an H-shaped pattern to 
avoid several very bright stars. The NTT Deep and 1511 fields contain tiles that are not 
contiguous. This was because the optical data were obtained with the sparse-array geometry 
Bernstein- Tyson Camera (BTC) on the 4-m telescope at CTIO. The journal of if-band 
imaging observations is given in Table 2, in which columns (1) through (5) and column (8) 
list the field, tile number, the 2000 coordinates of the tile centers, total exposure time, and 
date of observation, respectively. 

2.3. Optical Imaging Observations 

Optical images of the LCIR survey fields were obtained either from our own observations 
or from the literature using various large format CCDs at different telescopes. The journal of 
optical imaging observations is given in Table 3, which lists the field, telescope, instrument, 
plate scale, field of view (FOV), bandpasses, total exposure time, and date of observation. 
A more detailed description of the optical imaging observations will be presented elsewhere 
(Marzke et al. in preparation). 
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3. Image Processing 

In this section, we describe the procedures of processing the near-infrared //-band images 
of all six fields and summarize the qualities of the available optical images of the HDFS and 
CDFS regions. In addition, we discuss astrometric and photometric calibration of these 
images. 

3.1. The if-band Images 

The H-hand images obtained with CIRSI were processed using two independent re- 
duction packages, one developed by McCarthy, Wilson and Chen at Carnegie, the other by 
Lewis, Sabbey, McMahon, and Firth at the Institute of Astronomy (loA). 

The Carnegie reduction pipeline proceeds as follows. First, loop-combined images were 
formed by taking a mean of each loop of three to five exposures obtained without dithering 
the telescope. Aba rejection algorithm was applied to remove cosmic ray events, if there 
were more than four exposures in a loop. Next, the loop-combined images were divided by 
the normahzed dome flats and grouped together according to their locations on the sky. The 
use of a dome rather than dark-sky flat fleld image allows us to separate the pixel-to-pixel 
gain variations from the additive sky fringes which are non-negligible for these Hawaii arrays. 
The flat-field images have been normalized to remove the gain differences between the four 
detectors. These are generally less than 5%, except in the case of chip #3 for which the gain 
correction is ~ 20%. Next, a bad pixel mask was generated for each night using the flat- 
fleld exposures. Initial sky frames were constructed from medians of all frames for individual 
arrays that lay within a given pointing position. These sky frames were subtracted from each 
flat-fleld corrected image for the purpose of constructing an object mask for each image. The 
sky-subtracted images were registered using the offsets derived from the centroids of 5 to 15 
stars. All objects with peak intensities more than 10 a above the sky noise were identifled 
and sorted into four categories according to their isophotal sizes. Mask frames were then 
generated for each dithered exposure accordingly. Next, a sliding-median sky frame was 
obtained by taking the median of the preceding and succeeding three to flve images with 
the bad pixel and object masks applied. After being scaled by the ratio of the modes in 
the images, the median sky frame was subtracted from the appropriate flat-fleld-corrected 
image to remove the fringes due to OH emission lines from the sky that are present in the 
flat-fleld-corrected images. This latter scaling (typically ±1%) allowed us to remove non- 
hnear temporal variations in the sky level. Finally, the flat-fleld-corrected and sky-subtracted 
images were registered and combined using a 5 cr rejection criterion, a bad pixel mask, and a 
weight proportional to the inverse of the sky variance to form a flnal stacked image for each 
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pointing position. 

The contiguous 13' x 13' tiles were assembled from the 16 stacked images corresponding 
to the four arrays at four pointings of the telescope. The array boundaries were registered 
using stars in the overlap region, which was typically 20 to 30 arcsec in extent. The four 
arrays were not perfectly aligned with each other and a significant rotation (~ 0.4°) was found 
necessary to bring one of the arrays into registration with the other three. The final image 
was constructed by summing the 16 registered stacked images. A 5a rejection criterion and 
a weighting determined from the sky variance were also applied in the vicinity of the array 
boundaries. Ala error image was formed simultaneously for each tile through appropriate 
error propagation, assuming Poisson counting statistics. 

The loA based reduction pipeline is similar in conceptual design, although it is written 
in a somewhat different architecture and contains a few operational differences (Sabbey et 
al. 2001). The offsets between dithered images were derived via cross-correlations in the 
loA pipeline rather than from image centroiding. The object masks were derived from 
segmentation maps produced by SExtractor (Bertin & Arnout 1996) rather than using a set 
of five fixed object sizes as in the Carnegie pipehne. Finally, the assembly of the final tiles was 
completed after the application of a world coordinate system derived from comparison with 
the APM catalog. In the Carnegie pipeline the astrometry was derived for the completed 
tiles using the USNO astrometric catalog, discussed in § 3.3. 

The images of the first three tiles of the HDFS region, the first tile of the NTT field, and 
the second tile of the 1511 field contained significant artifacts associated with instabihties 
in the detector reset. These manifest themselves as low spatial frequency bias roU-offs that 
change shape and amplitude on short timescales. Both the Carnegie and loA pipelines deal 
with this low frequency signal by fitting low order polynomials to sigma-clipped images on 
a line-by-line and column-by-column basis. In the Carnegie pipeline this reset correction is 
applied after the final sky subtraction, in the loA pipeline it is applied before the fiat-field 
and sky corrections have been determined. New read-out electronics were implemented in 
August 2000, which largely suppress the reset instabihty. 

The full width at half maximum (FWHM) of a typical point spread function (PSF) and 
the 5 a detection threshold in a four arcsecond diameter aperture of each completed tile are 
given in columns (6) and (7) of Table 2, respectively. The PSF was measured to vary by 
no more than 10% across a tile, but the 5a sky noise was found to vary by as much as 50% 
across each tile. Figure 2 shows the distributions of the FWHM of a mean PSF (left) and 
the 5 a detection threshold estimated based on mean sky noise in a four arcsecond diameter 
aperture (right) for each completed tile. Figure 3 shows an example of a completed tile in 
the CDFS region from the Carnegie pipeline. The average PSF of this image has a FWHM 
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of 0.7 arcsec. The magnitudes of the two objects indicated by arrows are H — 19.4 and 
21.3. 



3.2. The Optical Images 

The optical U, B, V , R, and / images of an 40' x 40' area around the HDFS region were 
obtained and processed by a team at the Goddard Space Fhght Center during the HDFS 
campaign using the ETC (Tephtz et al. 2001). The spatial resolutions of the final combined 
images were measured to range from FWHM ?a 1.4 arcsec in the V band to ^ 1 . D arcsec 
in the B band, and the 5 a detection thresholds in a four arcsecond diameter aperture were 
measured to be [/ ~ 24.6, B ~ 26.1, V ~ 25.6, R ~ 25.1, and / ~ 24.2 magnitudes using 
the photometric zero points determined by the Goddard team. 

The optical V , R, I, and z' images of ~ 36' x 36' around the CDFS region were obtained 
using the Mosaicll imager on the 4-m telescope at CTIO in November of 1999. Individual 
images were processed following the standard mosaic image reduction procedures for data 
obtained with the Mosaicll imager, registered to a common origin, and coadded to form 
final combined images (Marzke et al. in preparation). The spatial resolutions of the final 
combined images were measured to range from FWHM ~ 1.0 arcsec in the / and z' bands 
to 1.2 arcsec in the V band, and the 5 a detection thresholds in a four arcsecond diameter 
aperture were measured to be F ~ 26.8, R ~ 26.2, I ~ 25.3, and z' ~ 23.7 mag. 

3.3. Astrometry 

Accurate astrometry is required for foUowup observations and for matching our cata- 
log with observations at other wavelengths. The optical images of the HDFS region were 
corrected for astrometric distortions with an RMS error of 0.16 pixels, corresponding to 
0.07 arcsec (Teplitz et al. 2001). We adopted the World Coordinate System (WCS) solution 
stored in the image headers to derive coordinates for objects identified in the images. To cor- 
rect astrometric distortions for the Mosaicll images of the CDFS region, we first re-sampled 
images of individual exposures to a Cartesian grid before co-addition and derived a best-fit, 
fourth-order polynomial for the WCS solution by matching stars in the U.S. Naval Observa- 
tory catalog (USNO-A2.0) with stars identified in the images (Marzke et al. in preparation). 
A total of 200 stars was included to obtain the best-fit astrometric solution. The goodness 
of fit may be characterized by an RMS residual of 0.2 arcsec. 

Astrometric solutions were also derived for the near-infrared images. We were able to 
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identify on average ~ 150 unblended and unsaturated stars in each if-band tile. Using a 
third order polynomial fit we were able to obtain solutions with rms residuals of between 
0.25 and 0.35 arcsec. 



4. Analysis 

In this section, we describe the procedures for object detection in individual bandpasses 
and object matching between frames taken through the various filters. We also examine the 
completeness limits of our detection algorithm and the reliability of photometric measure- 
ments in all bandpasses. 

4.1. Object Detection 

To obtain a uniform sample of distant galaxies of all types, it is crucial to identify 
galaxies of a wide range of optical and near-infrared colors. To reach this goal, we performed 
object detection using SExtractor for each bandpass in the regions where the if-band imaging 
survey was carried out. To reliably identify all the faint objects in each of the optical images, 
we set the SExtractor parameters by requiring that no detections be found in the negative 
image. First, we set the minimum area according to the FWHM of a mean PSF. Next, we 
adjusted the detection threshold to the lowest value that is consistent with zero detections 
in the negative image. To aid in deblending overlapping objects, we first constructed a 
"white hght" image by co-adding all the registered optical images that had been scaled to 
unit exposure time and filter throughput (see the next section for descriptions of image 
registration). Next, we applied SExtractor to the white light image and determined object 
extents based on their instrumental colors. Because the PSF varied between different optical 
bandpasses by at most 10% in each field, we were able to reliably deblend overlapping 
objects without erroneously splitting or blending individual objects in the white light image 
due to a large variation in the PSFs of different bandpasses. This procedure takes into 
account the intensity contrast between overlapping objects in all bandpasses instead of a 
single bandpass, allowing a more accurate object deblending for overlapping objects with 
different instrumental colors. 

We adopted a different approach to reliably identify all the faint objects in the if-band 
images. Because of a non-uniform noise pattern between individual arrays across a single tile, 
in particular in the join regions between adjacent array images, the zero negative detection 
criterion would force a significant underestimate of the survey depth in some areas. An 
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example to illustrate the noise pattern is given in Figure 4, which shows the associated 1 a 
noise image of the tile shown in Figure 3. We identified objects well into the noise in the H- 
band images and used our subsequent calibration of the completeness limits and photometric 
errors to set the faint magnitude limit of the i?-band catalog. 

Finally, to ensure accurate measurements of optical and near-infrared colors, we excluded 
objects that lie within 2 arcsec of the edges of each tile. 

4.2. Image Registration and Catalog Matching 

The optical U, B, V, R, and / images of the HDFS region were registered to a common 
origin by the Goddard team. The RMS dispersion in the object coordinates between different 
bandpasses was measured to be within 0.2 BTC pixels (corresponding to 0.08 arcsec). The 
optical V, R, I, and z' images of the CDFS region were registered individually for each area 
covered by a near-infrared tile, based on a solution determined using ~ 200 common stars 
identified in the tile. We applied standard routines to obtain the transformation solutions and 
to transform the images. We were able to register the optical images of the CDFS tiles by a 
simple transformation of low-order polynomials in the x and y directions without significantly 
degrading the image quality. To assess image degradation after the transformation, we 
measured the changes in the FWHM of the PSFs in transformed images and found an 
increase of no more than 5%. The RMS dispersion in the object coordinates between different 
bandpasses across an area spanned by a tile was measured to be within ~ 0.3 MosaicII pixels 
(corresponding to 0.08 arcsec). The dispersion was found to be the worst (as much as 0.5 
MosaicII pixels or 0.13 arcsec) in the join area between adjacent arrays of the MosaicII 
imager. To merge individual optical catalogs of each field, we stepped through each object 
identified in the white light image, examined the presence/absence of object detections at the 
location in individual bandpasses, and recorded the SExtractor measurements for detections 
and placed null values for no detections in the combined optical catalog. 

The H-band and optical catalogs were merged on the basis of astrometric position match- 
ing, rather than by transforming the H images onto the same pixels as the optical images. 
The large differences in pixel scales and in the PSFs of the optical and IR images result in 
significant image degradation if one transforms the H images to match the optical images. 
We derived a mapping between the H and optical coordinates from ~ 150 common stars in 
each tile. The RMS dispersion between the transformed optical coordinates and the if-band 
coordinates was measured to be ~ 0.2 CIRSI pixels (corresponding to 0.04 arcsec). Next, we 
repeated the catalog merging procedure described above by examining the presence/absence 
of object detections within a radius of two CIRSI pixels to the transformed object location 
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in the if-band image for all optically detected objects. The combined catalog was updated 
accordingly throughout the procedure. Finally, we compared objects in the combined optical 
catalog with those in the if-band catalog, identified those that appeared only in the i^-band 
catalog, and included them to form the final combined optical and near-infrared catalog for 
each tile. 



4.3. Survey Completeness 

Because detection efficiencies vary significantly for objects of different intrinsic profiles 
due to the strong surface brightness selection biases and because the noise patterns of the 
if-band images vary across individual tiles, particularly at the edges of individual arrays, 
it is crucial to understand the sensitivities and detection efficiencies of the LCIR survey in 
order to ensure the reliability of future analyses. To understand the completeness limits of 
the if-band images (similar analyses for the optical images will be presented in Marzke et al. 
in preparation), we performed a series of simulations, the results of which were also adopted 
for testing the photometric accuracy described in the next section. First, we generated test 
objects of different brightness for various PSF-convolved model surface brightness profiles. 
We adopted a Moffat (1969) profile to simulate the PSFs and considered (1) a point source 
model, (2) a de Vaucouleurs' r^l^ model characterized by a half-light radius r^, = 0.3 arcsec*, 
and (3) three exponential disk models characterized by scale radii of = 0.15, 0.3, and 0.6 
arcsec, respectively. The upper boundary of the selected disk scale lengths was determined 
according to the results of Schade ct al. (1996) based on galaxies identified in the CNOC 
survey. These authors found a mean scale length ~ 0.6 arcsec for field galaxies at 2; ~ 0.5. 
Next, we placed the test objects of a given magnitude at 1000 random positions in an H tile 
and repeated the object detection procedures using the same SExtractor parameters. Finally, 
we measured the recovery rate of the test objects as a function of H-hand magnitude. 

The results of the completeness test are presented in Figure 5 for Tile 1 of the CDFS field. 
Figure 5 shows that the 50% completeness hmit of the image is ~ 0.6 magnitudes fainter for 
an exponential disk of = 0.3 arcsec than for the largest disk model of rg — 0.6 arcsec and is 
an additional ~ 0.6 magnitudes fainter for a point source. The large difference between the 
completeness limits for objects of different intrinsic profiles makes an accurate incompleteness 
correction of the survey very challenging, because in practice we have no prior knowledge 
in the intrinsic profile of each source. It appears that under these circumstances correcting 
for the incompleteness based on one of these curves becomes arbitrary and therefore less 



^One arcsec corresponds to a projected proper distance of 5.6 h ^kpc at z = 1. 
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accurate. 

We repeated the simulations based on empirical surface brightness profiles determined 
from the image itself to address this problem and obtain a more meaningful estimate of 
the survey completeness. It is clear in Figure 5 that the completeness of an image is in 
practice set by the limit at which we can identify the majority of extended sources and 
that the survey for the largest disks in the image becomes incomplete at H at least 1.5 
magnitudes brighter than the 50% completeness limit for point sources. First, we adopted 
the stellarity index estimated by SExtractor, which according to Bertin & Arnouts (1996) 
is nearly unity for a point source and zero for an extended object. We identified objects of 
stellarity ^ 0.95 and with an if-band magnitude between H = 17 and 18, and coadded the 
individual images of these compact sources to form a high signal-to-noise empirical profile for 
point sources. Next, we repeated the simulations using the empirical profile to construct a 
completeness curve as a function of H for point sources. The results are shown as the thick, 
solid curve in Figure 5. Next, we determined the 50% completeness limit H^q for point 
sources and identified extended sources of stellarity ^ 0.5 and with an H-hand magnitude 
between H = H^q — 2 and H^^q — 1.5 (to include a large enough number of objects to form a 
smooth image). We co-added the individual images of these extended sources to form a high 
signal-to-noise empirical profile for extended sources. Finally, we repeated the simulations 
using the empirical profile to construct a completeness curve as a function of H for extended 
sources. The results are shown as the thick, dashed curve in Figure 5. The 95% completeness 
is ~ 0.5 magnitudes shallower for extended sources than for point sources. The completeness 
curve for the empirical profile of an extended source falls between those of Vg — 0.15 and 0.3 
arcsec, suggesting that sources with a large extended disk do not constitute a large portion 
of the faint population. We repeated this procedure for all H tiles. 



4.4. Photometric Accuracy 

Photometric cahbrations of the optical U, B, V, R, and / images were obtained by fit- 
ting a linear function, which included a zero-point, color, and extinction corrections, to the 
Landolt standard stars (Landolt 1992) observed on the same nights. Photometric calibra- 
tions of the z' images were obtained using the Sloan Digital Sky Survey secondary standards 
provided in Gunn et al. (2001), who adopted the AB magnitude system. We converted the 
best-fit zero-point to the Vega magnitude system, using z'p^^ = z' + 0.55 (Fukugita, Shi- 
masaku, & Ichikawa, 1995), for calibrating photometry of our objects. The RMS dispersion 
of the best-fit photometric solution was found to be ^ 0.03 magnitudes in all optical images 
of the CDFS field (a detailed discussion will be presented in Marzke et al. in preparation). 
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Photometric calibrations of the H-hernd images were obtained by fitting a constant term 
that accounts for the zero-point correction to a number of faint near-infrared standard stars 
(Persson et al. 1998) observed throughout every night. The RMS dispersion of the best-fit 
photometric solution was found to be ^ 0.02 mag. The zero-point was found to remain 
stable between different nights with a scatter of ^ 0.04 mag. 

To cross check the photometric accuracy of objects identified in the HDFS and CDFS 
regions, we first compared our photometry with that from the EIS project (da Costa ct al. 
1998) on an object-by- object basis. This clearly revealed an offset in the /-band photometry 
for the HDFS. Next, we compared optical colors of bright, compact sources (which arc likely 
to be stars) with the UBVRI colors of sub-dwarfs presented in Ryan (1989). We found that 
the stellar track agrees well with Ryan's measurements in all colors for objects in the CDFS 
region, but not for objects in the HDFS region in the U, V, and / bandpasses. Based on the 
results of stellar track comparison, we concluded that adjustments oi U — C/coddard + 0.3, 
V — Vcoddard + 0.1, and / = /coddard + 0.27 may be needed in order to bring the stellar 
colors into agreement, where f/coddard) K^oddardj and /ooddard ^-re the photometric zero-points 
published by the Coddard team at http://hires. gsfc. nasa.gov. /~ research/hdfs-btc/. This 
also brought our photometry very close to the EIS photometry of the same objects in the 
HDFS region. 

Accurate photometry of faint galaxies is hampered by the unknown intrinsic surface 
brightness profiles of the galaxies and because of strong surface brightness selection biases. 
Accurate total flux measurements of distant galaxies are crucial not only for an accurate 
assessment of galaxy evolution in luminosity function (e.g. Dalcanton 1998) but also for 
comparing the observed number counts with theoretical predictions based on different cos- 
mological models. Various approaches to measure the total fluxes of faint galaxies have been 
proposed. For example, an "adapted" aperture magnitude, also known as a Kron magni- 
tude, measures the integrated flux within a radius that corresponds to the flrst moment of 
the surface brightness distribution of an object and corrects for the missing flux outside the 
adopted aperture, which was found to be a constant fraction of the total flux for objects of 
different surface brightness proflles (Kron 1980). A simpler way is to measure the integrated 
flux within a flxed aperture and correct for the missing flux using a growth curve derived 
from the PSF of the image. 

To understand how well various approaches are able to recover the total flux for objects 
identifled in the survey, we performed a series of photometric tests for objects of different 
brightness and intrinsic surface brightness proflles based on the same simulations as described 
in § 4.3. First, we generated test objects of different brightness using a series of PSF- 
convolved model surface brightness profiles. Next, we placed the test objects at 1000 random 
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positions in an H tile and ran SExtractor on the simulated object image using the same 
photometric parameters. Finally, we compared various photometric measurements with the 
model fluxes. 

We present in Figure 6 comparisons of model fluxes and various photometric measure- 
ments for objects of different intrinsic surface brightness profiles to the 90% completeness 
limit of the image. The left most panel shows the test results for corrected isophotal mag- 
nitudes provided by SExtractor, the center panel for aperture magnitudes within a four 
arcsecond diameter, and the right panel for the "best" magnitudes estimated by SExtractor, 
which are either Kron magnitudes or corrected isophotal magnitudes in cases of crowded 
areas. The plotted points are the median residuals and the error bars indicate the 16th and 
84th percentiles of the residuals in 1000 realizations. Three interesting points are illustrated 
in Figure 6. First, the corrected isophotal magnitudes always underestimate the total flux at 
faint magnitudes, irrespective of object intrinsic surface brightness profile. Second, the four 
arcsecond diameter aperture photometry after a small, fixed aperture correction appears to 
deliver accurate total flux measurements for all but the largest disk model at all magnitudes 
to the 90% completeness limits. Third, the SExtractor "best" magnitudes appear to be able 
to recover the total flux of every profile except for the largest disks at the faintest magni- 
tudes, although it also appears that a small amount of "aperture" correction is needed to 
bring the centroids of the residuals to zero. Because galaxies of scale length greater than 0.5 
arcsec, corresponding to 0.8 arcsec half-hght radius, appear to be rare (e.g. Yan et al. 1998), 
we conclude that the corrected four arcsecond aperture magnitude provides the most accu- 
rate estimate of the total object flux. Although the use of a relatively large aperture results 
in reduced signal-to-noise ratios and hence a loss in depth, the loss in precision is offset by 
the improved accuracy in both the total magnitude measurements and the determination of 
colors from images with a wide range of seeing. 



4.5. Star and Galaxy Separation 

Stellar contamination of the final galaxy catalog from the LCIR survey is non-neghgible, 
particularly in the HDFS. An accurate separation of stars from the galaxy catalog is therefore 
essential for subsequent statistical analysis. To identify stars in the survey fields, we can 
either rely on morphological classification by comparing object sizes and shapes with the 
PSF, or we can classify objects on the basis of their spectral shapes. Neither approach by 
itself can robustly identify stars at all signal-to-noise levels. 

To test the reliability of the stcllarity index provided in SExtractor at different signal 
levels, we measured the rate at which a point source has a measured stellarity index less 
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than a threshold value in the completeness test described in § 4.3. The results are shown 
in Figure 7. The solid curve indicates the fraction of the input stars that have a stcllarity 
index <^ 0.5 as a function of //-band magnitude. The dashed curve indicates the fraction 
of stars that have a stellarity index <^ 0.8. The dotted line indicates the 50% completeness 
limit (H^q) of the test image for a point source. We found that while SExtractor is able to 
uncover point sources at fairly bright magnitudes, it becomes unreliable at around H^q — 2 
magnitudes (i.e. more than 20% of stars would be misclassified as galaxies) and unstable 
at around H50 — 1.5 magnitudes (where objects are more or less equally divided between 
compact and extended sources irrespective of their true shapes), and therefore cannot by 
themselves be adopted to robustly separate stars from galaxies. In addition, identifying 
stars based on morphological classification alone introduces biases against compact sources 
such as quasars, narrow emission line galaxies, etc. 

To test the reliability of comparing object colors with known stellar energy distributions, 
we first included a suite of empirical stellar templates collected from the literature (Oppen- 
heimer et al. 1998; Pickles 1998; Fan et al. 2000; Lcggctt ct al. 2000) in the photometric 
redshift likelihood analysis (Chen et al. in preparation). The adopted stellar templates span 
a wide range in both spectral coverage, from optical to near-infrared wavelengths, and stel- 
lar types, from early-type O and B stars to late-type L and T dwarfs. Then we compared 
the results of photometric redshift analysis with spectroscopic redshift measurements col- 
lected from the literature for objects in the HDFS region. We found that while all the 
spectroscopically identified stars were accurately identified as stars based on the photomet- 
ric redshift analysis, ~ 10% of the spectroscopically identified galaxies were misidentified as 
stars. Therefore, we conclude that photometric redshift analysis with our present filter set 
alone cannot robustly separate galaxies from stars. 

To improve the accuracy of star and galaxy separation, we investigated a method of 
combining both approaches. We first examined the stellarity index estimated for each object 
by SExtractor in all bandpasses. For objects with a stellarity index greater than 0.95 in one 
or more bandpasses, we used both stellar and galaxy templates in the photometric redshift 
analysis. For objects with a stellarity index less than 0.95 in all bandpasses, we used only 
galaxy templates in the photometric redshift analysis. Comparison of photometric and spec- 
troscopic redshifts for objects in the HDFS region demonstrates that all 11 spectroscopically 
identified stars may be reliably identified in this procedure. 

The addition of other colors (e.g. J & K) will improve the photometric separation of 
stars and galaxies. At present, the star and galaxy separation may be inaccurate at faint 
magnitudes where the stellarity index becomes unreliable. However, galaxies are likely to 
dominate the number counts at faint magnitudes (see § 6.2). Figure 8 shows the stellar frac- 
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tion in the HDFS and CDFS regions based on the results of the star and galaxy separation. 
Closed circles indicate the mean surface density of galaxies observed in the two fields (see 
§ 6.2 for a detailed discussion); stars indicate the surface density of stellar objects in the 
HDFS region; and open triangles indicate the surface density of stellar objects in the CDFS 
region. It shows that the star counts increase gradually with the H magnitude, while the 
galaxy counts show a much steeper slope. The smooth, shallower slope of the star counts 
suggests that stars do not contribute a significant fraction at if ;^ 19 in both HDFS and 
CDFS fields. 



5. The Catalogs 

The if-band imaging survey presently covers 1.1 square degrees of sky to a mean 5 a 
detection limit in a four arcsecond diameter aperture of if ~ 20.8. Here we present sample 
near-infrared and optical photometric measurements for objects identified in the HDFS and 
CDFS regions. We have performed object detection in all bandpasses and identified ~ 24, 000 
galaxies over 847 arcmin^ in the HDFS region, of which ~ 6, 700 are detected in the if-band 
survey, and ~ 30, 000 galaxies over 561 arcmin^ in the CDFS region, of which ~ 7, 400 are 
detected in the if-band survey. Among the if-detected galaxies in the two regions, ~ 2, 000 
have I — H > 3 and are likely to be primarily evolved galaxies at redshifts z > 1. 

The catalogs for each field is organized as follows: For each object we give an identifica- 
tion number, x and y pixel coordinates in the optical images, a (J2000) & () (J2000), the Kron 
radius (in arcseconds) for each filter, the stellarity index determined by SExtractor, 2 " and 
4" diameter aperture magnitudes and their associated uncertainties, the corrected isophotal 
magnitude, auto magnitude and the best magnitude and their assocaited uncertainties as pro- 
vided by SExtractor for all the objects detected in each filter. We present an example of the 
object catalog in Table 4 and an example of the photometric catalog in Table 5, which list the 
first 50 objects identified in the CDFS region. Complete near-infrared and optical photomet- 
ric catalogs of the two fields may be accessed at http://www.ociw.edu/lcirs/catalogs.html. 
Catalogs for the remaining fields and their completeness limits will be posted at this address 
as they become available. 

6. Discussion 

As we identify objects in each optical and near-infrared bandpass the catalogs contain 
objects with a wide range of colors, ranging from the reddest objects, whether they be old 
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or dusty, to the most actively star forming systems. The full galaxy sample is therefore 
representative of galaxies of a wide range of stellar properties and redshifts, and may be 
used to address a variety of cosmological issues. Here we examine photometric properties 
of the if-band selected objects and derive the number-magnitude relation for the if-band 
selected galaxies as well as for red galaxies of different optical and near-infrared colors. We 
compare our measurements with various model predictions and discuss the implications. 

6.1. Distributions of Optical and Near-infrared Colors 

We present the I — H versus if-band color-magnitude diagram for the ~ 14, 000 i7-band 
selected galaxies identified in the HDFS and CDFS regions in Figure 9. Objects classified 
as stars are not included. We show the 95% completeness limit of the if-band data in the 
short dashed fine and the 5 a detection threshold over a four arcsecond diameter aperture for 
the /-band images of the HDFS and CDFS regions in the dot-dashed fines. We also present 
the expected color-magnitude evolution of an L* elliptical galaxy from 2; = 2 to the present 
time under different evolution scenarios for h — 0.7. The curves trace the trajectories of 
a non-evolving elliptical galaxy spectrum (blue), a passively evolving galaxy formed in a 
single burst at 2; = 10 (red), and a galaxy formed at 2; = 30 with an exponentially declining 
star formation rate (SFR ~ exp(— t/r)) for r = 1 Gyr (green). The blue dashed curve 
corresponds to a non-evolving 3 early-type galaxy. 

Comparisons of measurements and model predictions show that the I — H colors of 
galaxies with ^ 18 agree well with the expected colors at low redshifts. In addition, they 
show that the median color becomes progressively bluer at fainter magnitudes, accompanied 
by a tail of galaxies with very red colors {I — H ^ 3). These red galaxies constitute a 
negligible fraction of the total population at if ~ 18, but become a significant constituent 
by if ~ 20. Similar trends are seen in B — K, R — K and I — K color-magnitude diagrams 
based on smaller area surveys (Elston et al. 1988; Cowie et al. 1994; EUis 1997; Cowie et al. 
1997, Menanteau et al. 1999). Furthermore, the good agreements between observed colors 
and model predictions suggest that these red galaxies are mildly evolving early-type galaxies 
at 2; ~ 1 sampled over the bright end of the luminosity function (~ 1 — 5L*). The red 
sequence observed in clusters with z ~ 0.8 — 1.2 populate a similar range in the I — H vs. H 
color-magnitude diagram, supporting the idea that the bulk of the red galaxies are evolved 
early-type galaxies at intermediate redshifts (e.g. Stanford et al. 1997, 1998; Chapman, 
McCarthy & Persson 2000). 

We present the V ~I vs. I—H color-color distribution of the i/-band detected objects in 
the HDFS and CDFS regions in Figure 10. Stars that were classified according to the criteria 
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described in § 4.5, are indicated by crosses and galaxies are indicated by filled circles. Open 
circles indicate galaxies that are not detected in V. The solid curves trace the trajectories of 
different CWW spectral templates at redshifts from z = 2 to the present time: elliptical/SO 
(yellow), Sab (green), Scd (cyan), and irregular (blue). The dotted curve shows the trajectory 
of a passively evolving galaxy formed in a single burst at redshift Zf = 10. The short and 
long dashed curves trace the colors of galaxies formed aX zj = 30 with exponentially declining 
star formation rates for r = 1 and 2 Gyr, respectively. The dash-dotted hne corresponds to 
a galaxy formed at Zf — 30 with a constant star formation rate. Filled circles along each 
curve indicate the corresponding redshifts, starting from z — 2 in the lower right corner, 
evolving toward lower redshifts in steps of Az = 0.5, and ending at 2; = in the lower left 
corner. 

The general agreement between distributions of various model predictions and measure- 
ments of galaxies in our catalog, and the clear separation of the stellar sequence suggest that 
most stars have been accurately identified in these two survey fields. A similar separation 
is made in a i? — / vs. I — K color-color diagram (e.g. Gardner 1992; Huang et al. 1997). 
The large scatter in the V — I colors spanned by the galaxies reveals a wide range of star 
formation histories for which our simple models provide only an idealized description. In 
particular, the large scatter in the V — I colors for galaxies of I — H colors indicative of 
redshifts 2; ~ 1 suggest a significant amount of ongoing star formation in these galaxies. 
Similar results are found in detailed studies of morphologically selected early-type galaxies 
at intermediate redshifts that reveal a deficit of distant ellipticals with optical and UV colors 
that match with expectations based on passively evolving galaxies formed at high redshifts 
(e.g. Menantcau ct al. 1999; Menanteau, Abraham, & Ellis 2001). In addition, the small 
number (< 15%) of galaxies with I ~ H ^ 3 and V — I ^ 3 suggests that there are few 
objects that may be adequately characterized as pure passively evolving systems formed at 
high redshifts. The trajectories of the evolving galaxy models in the color-color diagram 
suggests that these red galaxies span a significant range in redshift and that those of bluer 
V — I colors largely arise at redshifts z ^ 1.5. This is supported by their angular clustering 
amphtude (McCarthy et al. 2001b) and photometric redshift measurements (Chen et al. in 
preparation) . 



6.2. The i7-band Number Counts 

Galaxy number counts provide a diagnostic for distinguishing between different evolu- 
tionary scenarios for galaxies of different colors (e.g. Metcalfe et al. 1996). Comparisons of 
number count measurements from different surveys also allow us to examine the accuracy of 
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our photometric calibration. Figure 11 shows the if-band differential galaxy counts N{in) 
derived from ~ 14, 000 galaxies selected in the if-band survey toward the HDFS and CDFS 
regions (filled circles). We applied incompleteness corrections down to the 95% completeness 
level using the incompleteness curves determined for compact and extended sources in each 
tile. According to the procedures described in § 4.3, we adopted the stellarity index deter- 
mined in SExtractor for each object and corrected for the incompleteness at faint magnitudes 
using the compact source incompleteness curve for objects of stellarity index greater than 
0.5 and using the extended source incompleteness curve otherwise. The vertical bars indi- 
cate the uncertainties in our measurements estimated using a variation of the "bootstrap" 
resampling technique. The bootstrap resampling technique explicitly accounts for the sam- 
pling error and for photometric uncertainties. Wc first rcsamplcd objects from the parent 
catalog (allowing for the possibility of duplication) and added random noise (within the 
photometric uncertainties) to the photometry. Then we estimated the H-hand differential 
number counts for the resampled, perturbed catalog. Finally, wc repeated the procedures a 
thousand times and determined the 1 a uncertainties based on the 16th and 84th percentiles 
of the distributions in each magnitude bin. 

For comparison, we have also included in Figure 11 previous iJ-band measurements 
reported by various groups. Our measurements span six magnitudes in H and nicely bridge 
the gap between the existing measurements at the bright (open squares from Martini 2001b) ^ 
and faint ends (Yan et al. 1998). The excellent agreements in both the slopes and the normal- 
ization demonstrate not only the accuracy of our iJ-band photometry, but also the accuracy 
in our star and galaxy separation. The best fitting power-law slopes, d\og N{m.)/dm, were 
found to be 0.45 ± 0.01 mag-^ at ^ 19 and 0.27 ± 0.01 mag"^ at H ;^ 19 with a mean 
surface density of ~ 7, 200 dcgrcc^^ mag^^ at H = 19^. The slope of the faint-end iJ-band 
counts is in line with that of the faint-end X-band counts (Cowic et al. 1994; Djorgovski et 
al. 1995; Bcrshady, Lowcnthal, & Koo 1998; Metcalfe ct al. 1996; McCracken et al. 2000). 
The slope of the bright-end if-band counts {H <^ 19) is slightly shallower than that of 
the bright-end X-band counts reported by Huang et al. (1997) and Metcalfe et al. (1996), 
although we note that our counts are not well determined at H < 16. 



^Part of the apparent excess in the number count measurements by Martini (2001b) is due to an additional 
aperture correction that the author apphed in the galaxy photometry to account for missing Hght in the 
wings of galaxy profiles, which was found to be between 0.1 and 0.2 mag. 

^We note that the uncertainties in the number count measurements due to random errors are small 
given the large sample from our survey, but sysmatic errors are likely to dominate and more difficult to 
measure. In particular, number count measurements at the faintest magnitude bins are more likely to be 
overestimated due to photometric undertainties of the even fainter and more abundant populations. As a 
result, the faint-end slope is likely to be overestimated. 
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To compare with simple model predictions, we considered a no-evolution model based 
on an early- type galaxy template (the solid curve), a passive evolution scenario for galaxies 
formed in a single burst at Zf = 10 (the dotted curve), an exponentially declining star 
formation for galaxies formed at Zf = 30 (the short dashed curve), and a constant star 
formation for galaxies formed at Zf = 30 (the long dashed curve). We adopted the K- 
band luminosity function derived by Gardner et al. (1997) and applied a.n H — K ^ 0.22 
color correction (Mobasher, Ellis, & Sharpies 1986) to obtain an estimate for Mh^. This 
is appropriate for galaxies of different type because the /c-correction is nearly independent 
of spectral type at near-infrared wavelengths. We find that all but the passive evolution 
scenario (the dotted curve) agree well with the measurements. 



6.3. Surface Density of Red Galaxies 

In this section, we present surface density estimates of moderately red galaxies of I—H ^ 
3 and compare the measurements with those of extremely red objects (EROs) identified in 
the survey. EROs are typically selected based on their extreme optical and near-infrared 
colors i? - X ;^ 6 (e.g. Elston, Rieke, & Rieke 1988; McCarthy, Persson, & West 1992; Hu & 
Ridgway 1994). Here we selected EROs using R — H ^ 5, which according to the simulations 
presented in Figure 1 is equivalent to R — K ^6 for high-redshift evolved populations. 

Figure 12 shows the if-band differential number counts for galaxies of / — if ;^ 3 (filled 
circles) and for galaxies oi R — H ^ 5 (squares) identified in the HDFS and CDFS regions. 
Incompleteness corrections were applied at faint magnitudes and the uncertainties indicated 
by the vertical bars were estimated following the same procedures described in § 6.2. The 
slope of the red counts is quite steep, rflog N{m)/dm = 0.84 ± 0.06 mag~^ for galaxies of 
I — H ^3 and 0.89 ± 0.24 mag~^ for galaxies of — if ;^ 5 at if ^ 20. The counts appear 
to have shallower faint-end slopes for both of the red subsamples — between 0.17 and 0.32 
mag~^ — at H ^ 20. While both subsamples show consistent slopes at the bright end, they 
have a more steeply rising N{m) with the H magnitude than the total if-band detected 
population. The steeper bright-end slope of the red galaxies indicates that most of the 
foreground {z <^ 1) galaxies have been effectively excluded from the red samples and that 
the shape of the number counts reflects the shape of the bright-end luminosity function of the 
red population. In addition, the consistent slopes between the two red subsamples suggest 
that the underlying luminosity function is very similar. 

For comparison, we also calculated the predicted differential number counts for red 
galaxies with I — H ^3 based on the no-evolution scenario and the exponentially declining 
star formation scenario with r = 1 Gyr that best fit the observed total if -band number 
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counts, as discussed in § 6.2. To model the intrinsic luminosity function of evolved galaxies, 
we adopted the population ratio of elliptical galaxies in the local iinivcrsc (e.g. Binggeli, 
Sandage, & Tammann 1988). We first scaled the i^T-band luminosity function obtained by 
Gardner et al. (1997) for the total galaxy population by M^{red) = M^{tot) — 0.2 and 
(f)*{red) = 0.2 (f)^,{tot) and then assumed a faint-end slope of a = 1. The curves indicate 
the results of the calculations and are the same as those in Figure 11. It appears that both 
scenarios fit the observations fairly well, but that the exponentially declining star formation 
model predicts a sharp turnover in the differential number counts at H ^21. The turnover 
is due to the fact that at these faint magnitudes the evolved population is dominated by 
those at redshifts 2; ;^ 2 where they would have bluer I — H colors under the exponentially 
declining star formation scenario and would therefore be excluded based on the I ~ H ^ 3 
selection criterion (sec Figure 9). Because of large uncertainties in the optical and //-band 
photometry for galaxies at these faint magnitudes, we cannot rule out either of the scenarios. 
The moderately good agreements between observations and model predictions, together with 
the large scatter in the optical V — I color discussed in § 6.1, however suggest that there is 
still some amount of star formation going on in the red population at redshifts z > 1. 

In Figure 13 we plot the cumulative surface density measurements of galaxies with 
I — H ^3 (filled circles) and of R — H ^ 5 (squares) as a function of if-band magnitude, in 
comparison with those of the total i^-band selected galaxies and with previous measurements 
of EROs presented by Daddi et al. (2000; green squares), Yan et al. (2000; blue triangles), 
and Thompson et al. (1999; red cross). Our surface density measurements of EROs span 
two magnitudes in H and are consistent with those from previous surveys. In addition, we 
find that red galaxies with I — H ^ 3 constitute ~ 20% of the iJ-band detected galaxies 
at H <^ 21, but only ^ 2% at H <^ 19. Specifically, we expect to observe on average ~ 40 
galaxies of I — H ^3 and ~ 8 galaxies of R — H ^5 among ~ 270 iJ-band detected galaxies 
at ^ 21, but only ~ 1 galaxy of I — H ^3 and less than one of R — H ^ 5 among ~ 46 
if-band detected galaxies at iJ ^ 19 in a random 25-arcmin^ region. 

On the other hand, because EROs and these moderately red galaxies are found to exhibit 
strong angular clustering (Daddi et al. 2000; McCarthy et al. 2001b; Firth et al. 2001), a 
large disperson in the surface density measurements is expected. To examine the accuracy 
of the surface density measurements presented in Figure 13, we estimated the field to field 
variation by measuring the surface density of galaxies of / — if ;^ 3 in a randomly selected 
5x5 arcmin^ region in each tile and compared the measurements between different tiles. 
We found that the total number of red galaxies of H ^ 20.5 varies from 15 to 40 with 
a mean surface density of ~ 1 arcmin"^ and an rms dispersion of ~ 0.4 arcmin"^. The 
mean surface density estimated based on the randomly selected smaller areas is consistent 
with the measurements presented in Figure 12 and the red counts over 0.62 square degrees 
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presented in McCarthy et al. (2001a,b), indicating that our surface density measurements 
are not significantly affected by the cosmic variance due to clustering. In addition, following 
Roche et al. (1999), we estimate an angular clustering amplitude of these red galaxies to be 
A ^ A aX 1" (corresponding to a correlation angle of ~ 5 arcsec), which indeed agrees 
well with the results of an angular clustering analysis presented in McCarthy et al. (2001b). 
We therefore conclude that a deep, wide-field survey is necessary not only to identify a large 
number of red galaxies, but also to obtain accurate estimates of their statistical properties. 

7. Summciry 

We present initial results from the Las Campanas Infrared Survey, which is a wide- 
field near-infrared and optical imaging survey designed to reliably identify a statistically 
significant and representative sample of evolved galaxies at redshifts z > 1 and to secure a 
uniform sample of galaxies of all types at intermediate redshifts. We present sample near- 
infrared and optical photometric catalogs for objects identified in two of the six survey fields, 
surrounding the Hubble Deep Field South and Chandra Deep Field South regions. We have 
identified ~ 14,000 i?-band detected galaxies over 1,408 arcmin^ of sky, ~ 2,000 of which 
are likely to be evolved galaxies at 2; ;^ 1 with I — H ^3. Furthermore, we find that: 

1. Red galaxies with I — H ^ 3 are a significant constituent of galaxies with H ^ 20 
and are consistent with mildly evolving early- type galaxies at 2; ~ 1. But there exists a large 
scatter in the V — I colors spanned by the red galaxies, indicating a significant amount of 
on-going star formation in these galaxies. 

2. The differential number counts N(m) for the i?-band detected objects has a slope of 
dlog N{m)/dm = 0.45 ± 0.01 mag'^ at F 19 and 0.27 ± 0.01 mag-^ at ;^ 19 with a 
mean surface density ^ 7, 200 degree"^ mag~^ H — 19. 

3. The differential number counts for red objects (/ — if ;^ 3) detected in the if-band 
has a very steep slope, dlog N{m;I — H ^ 3) /dm — 0.84 ± 0.06 mag~^ at if ^ 20 and 
0.32 ± 0.07 mag~^ a,t H ^ 20, with a mean surface density fa 3, 000 degree"^ mag~^ at 

4. Red galaxies {I — H ^ 3) constitute 20% of the if-band detected galaxies at 
H ^ 21, but only 2% at H ^ 19. This suggests that foreground galaxies {z ^ 1) may be 
effectively excluded from the total i7-band population by apply the I — H color selection 
criterion and that the steep bright-end slope of red galaxy number counts reflects the shape 
of the bright-end luminosity function of the red population. 
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5. Red galaxies are strongly clustered, producing strong field to field variation in their 
surface density. We estimate a mean surface density of ^ 1 arcmin^^ with an rms dispersion 
of ~ 0.4 arcmin"^ over 5x5 arcmin^ sky regions for galaxies with H <^ 20.5 and I — H ^ 3, 
implying an angular clustering amplitude of ^4 fti 4 at 1" and a correlation angle of ~ 5 
arcsec. 

We are grateful to Xiaohui Fan, Sandy Leggett, and Ben Oppenheimer for providing 
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Campanas Observatory. This research was supported by the National Science Foundation 
under grant AST-9900806. The CIRSI camera was made possible by the generous support 
of the Raymond and Beverly Sackler Foundation. 
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Fig. 1. — The predicted H-hand magnitude (left panel) and optical and near- infrared colors 
(right panel) of an L* elliptical galaxy as a function of redshift. The solid curves are for the 
non-evolving model; the dashed curves are for a single burst galaxy formed at Zf = 10; the 
dotted curves are for a single burst galaxy formed at Zf = 5. The curves of the right panel 
represent the redshift evolution of (from top to bottom) the R — K, R H, I — H, and 
z' — H colors. It is shown that evolved galaxies at rcdshifts z > 1 may be selected based on 
R-K ;>6, R-H ^ 4.5, I - H ^ 3, or z' - H ;> 2. 

Fig. 2. — Summaries of the image quality of the if-band survey. The left panel shows the 
distribution of the FWHM in arcsec of a mean PSF in each tile. The right panel shows the 
distribution of the 5 a H-hand detection threshold in each tile estimated based on mean sky 
noise over a four arcsecond diameter aperture. 

Fig. 3. — An example of an assembled tile in the CDFS region. The image is approximately 
13.5 arcmin on a side. The 5 a point source detection limit is iJ ~ 20.9 in a four arcsecond 
diameter aperture in this tile and the PSF is ~ 0.7 arcsec FWHM. The bottom panel shows 
details of an 2 x 1 arcmin^ region. To illustrate the depth and contrast of the image, we 
have selected two objects oi H = 19.4 (left) and 21.3 (right), respectively, as indicated by 
the arrows. 

Fig. 4. — The associated 1 a noise image of Figure 3. The image dimension is the same 
as Figure 3. Darker regions represent noiscr areas. The composition of 16 stacked images, 
corresponding to four arrays at four pointings of the telescope, is clearly illustrated by the 
"square" pattern in the image. The noise level varies between arrays and pointings because of 
the differences in array sensitivity and total exposure time. Black dots indicate the locations 
of bright stars, where the noise level is dominated by object flux instead of sky. 

Fig. 5. — The detection efficiency estimated in one of the CDFS H tiles for objects of 
different intrinsic profiles and angular sizes. As indicated in the upper-left corner, the solid 
curve represents a stellar profile; the dotted curve represents an r^^^ de Vaucoulcur profile 
with a half-hght radius of 0.3 arcsec; the long-dashed, short-dashed, and dash-dotted curves 
represent disk profiles with scale radii of 0.15, 0.3 and 0.6 arcsec, respectively. The thick 
curves show the detection efficiencies for a point source profile (solid) and an extended profile 
(short-dashed) determined empirically by coadding individual object images (see § 4.3 for a 
detailed description). 

Fig. 6. — Comparisons of model fluxes and various photometric measurements for objects 
of different brightness and intrinsic surface brightness proflles. The flrst column shows the 
test results for corrected isophotal magnitudes provided by SExtractor, the second column 
for aperture magnitudes within a four arcsecond diameter, and the third column for the 
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"best" magnitudes estimated by SExtractor. Points are the median residuals and error bars 
indicate the 16th and 84th percentile of the residuals in 1000 realizations. The last point in 
each panel corresponds to the 90% completeness limit of the image, which varies with object 
intrinsic surface brightness profile. 

Fig. 7. — Robustness of using the stellarity index provided in SExtractor to identify stars. 
According to Bertin & Arnouts (1996), a point source has a stellarity index close to one, while 
an extended object has a stellarity index close to zero. The solid curve indicates the fraction 
of the input stars that have a stellarity index ^ 0.5 as a function of H-hand magnitude. The 
dashed curve indicates the fraction of stars that have a stellarity index <^ 0.8. The dotted 
line indicates the 50% completeness limit of the test image for a point source. 

Fig. 8. — Stellar fraction in the HDFS and CDFS regions. Closed circles indicate the mean 

surface density of galaxies observed in the two fields; stars indicate the surface density of 
stars in the HDFS region; open triangles indicate the surface density of stars in the CDFS 
region. 

Fig. 9. — The I—H color versus i?-band magnitude diagram for the if-band selected galaxies 
in the HDFS and CDFS regions. We show the 95% completeness limit of the if-band data in 

the short dashed line and the 5 a detection threshold over a four arcsecond diameter aperture 
for the /-band images of the HDFS and CDFS regions in the dot-dashed lines. The blue 
solid curve represents a non-evolving elliptical galaxy; the red curve represents a passively 
evolving galaxy formed at z = 10; and the green curve represents an L^, galaxy constantly 
forming stars since z = 30. The dashed blue curve represents a non-evolving 3L^, elliptical 
galaxy. The dotted hue indicates the I — H — 3 color selection threshold for red galaxies at 

Fig. 10. — The V — I vs. I — H color-color diagram for objects identified as stars (crosses) 
and galaxies (filled points) in the HDFS and CDFS regions. Open circles indicate galaxies 
that are not detected in V. Solid curves indicate the mean colors of galaxy of different types 
as a function of redshift at redshifts between z = to 2, from elliptical or SO (yellow). Sab 
(green), Scd (cyan), to irregular (blue) galaxies. The dotted curve indicates the predicted 
colors of a passively evolving galaxy formed in a single burst at redshift Zf — 10. The short 
and long dashed curves indicate the predicted colors of galaxies formed at Zf — 30 with an 
exponentially declining star formation rate exp(— t/r) for r = 1 and 2 Gyr, respectively. The 
dash-dotted line indicates the predicted colors of a galaxy formed at 2;/ = 30 with a constant 
star formation rate at all redshifts. Filled circles along each curve indicate the corresponding 
redshifts, starting from z = 2 on the right and evolving toward lower redshifts in steps of 
Az — 0.5. The dotted line indicates the I — H — 3 color selection threshold for red galaxies 
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at z ;^ 1. 

Fig. 11. — Differential number counts for the if-band selected galaxies in the HDFS and 
CDFS regions. The open symbles are previous measurements presented by Yan et al. (1998; 
triangles) and Martini (2001b; squares). The curves are predictions based on various star 
formation scenarios discussed in § 6.2. 

Fig. 12. — Differential number counts for the H-hand selected red galaxies in the HDFS 
and CDFS regions. Points are for galaxies oi I — H ^ 3 and squares are for galaxies of 
R — H ^ 5. As shown in Figure 11, the solid curve indicates the prediction based on a 
non-evolving scenario and the short-dashed curve indicates an exponentially declining star 
formation scenario for galaxies formed at Zf — 30. These models are described in § 6.3. 

Fig. 13. — Cumulative surface densities of the total if -band detected galaxies (open circles) 
and galaxies of I — H ^ 3 (points) and of R — H ^ 5 (squares) as a function of if-band 
magnitude. The green squares are measurements presented by Daddi et al. (2000); the blue 
triangles are from Yan et al. (2000); and the red cross is from Thompson et al. (1999; red 
cross) . The solid hues are the best-fit power laws to the cumulative counts. 
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Table 1. LCIR Survey Fields 



Field 


RA (2000) 


DEC (2000) 




NOAO-Cetus 


02^ 


10™ 


-04° 


37' 


60° 


CDFS 


03^ 


32™ 


-27° 


38' 


54° 


NTT Deep . . . 


12*^ 


05™ 


-07° 


27' 


52° 


1511 


15'^ 


24™ 


+00° 


10' 


44° 


SSA22 


T2^ 


00™ 


+00° 


21' 


40° 


HDFS 


22^ 


33™ 


-60° 


39' 


49° 
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Table 2. Journal of i?-band Imaging Observations 



Exposure FWHM 



Field 


Tile 


RA (2000) 


DEC (2000) 


(sec) 


(arcsec) 


H5a^ 




Date 




(1) 


(2) 


(3) 


(1) 


(■^) 


(6) 


(7) 











NOAO-Cetus 


1 


02:10:13.0 


-04:33:00.0 


4800 


0.80 


20.1 


Oct 


16- 


-19 


1999 




2 


02:09:22.0 


-04:33:00.0 


2430 


0.70 


20.2 


Oct 


09-09 


2000 




3 


02:10:13.0 


-04:45:48.0 


3300 


0.72 


20.2 


Oct 


11- 


-12 


2000 




4 


02:09:22.0 


-04:45:48.0 


2320 


0.68 


20.3 


Oct 


09- 


-09 


2000 


CDFS 


1 


03:32:36.0 


-27:46:30.0 


4420 


0.66 


21.3 


Oct 


07- 


-09 


2000 




2 


03:32:36.0 


-27:33:53.0 


4860 


0.70 


21.2 


Oct 


05- 


-07 


2000 




3 


03:31:39.0 


-27:33:53.0 


5535 


0.69 


21.2 


Sep 


15- 


-16 


2000 




4 


03:31:39.0 


-27:46:30.0 


6243 


0.80 


20.8 


Dec 


20- 


-31 


1999 


NTT Deep . . . 


1 


12:05:57.8 


-07:41:58.0 


4700 


0.80 


20.5 


Feb 


13- 


-16 


1999 




2 


12:05:57.8 


-07:21:52.0 


4860 


0.90 


20.9 


Jan 


01- 


-17 


2000 




3 


12:04:36.6 


-07:41:58.0 


4691 


0.80 


20.6 


Feb 


13- 


-16 


2000 




4 


12:04:36.6 


-07:21:52.0 


5264 


0.75 


21.3 


Dec 


26- 


-30 


1999 


1511 


2 


15:24:07.7 


+00:20:16.0 


3600 


0.68 


20.1 


Feb 


13- 


-16 


2000 


SA22 


1 


22:18:00.0 


+00:14:34.0 


4321 


0.66 


20.1 


Oct 


16- 


-24 


1999 




2 


22:17:08.0 


+00:14:34.0 


3510 


0.60 


20.8 


Oct 


06-07 


2000 




3 


22:18:00.0 


+00:27:30.0 


4620 


0.77 


20.8 


Oct 


04- 


-06 


2000 




4 


22:17:08.0 


+00:27:30.0 


3731 


0.65 


20.5 


Oct 


12- 


-13 


2000 


HDFS 


1 


22:33:13.0 


-60:39:27.0 


4826 


0.92 


20.5 


Oct 


18- 


-19 


1999 




2 


22:34:56.0 


-60:39:27.0 


3645 


0.90 


20.5 


Oct 


20- 


-29 


1999 




3 


22:31:30.0 


-60:39:27.0 


4758 


1.06 


20.6 


Oct 


25- 


-29 


1999 




4 


22:34:56.0 


-60:26:50.0 


3645 


0.84 


21.3 


Sep 


07-07 


2000 




5 


22:34:56.0 


-60:49:50.0 


4230 


1.00 


21.1 


Sep 


13- 


-15 


2000 




6 


22:31:30.0 


-60:49:50.0 


6986 


0.85 


21.2 


Oct 


04- 


-11 


2000 




7 


22:31:30.0 


-60:26:00.0 


4890 


0.85 


21.3 


Oct 10-12 2000 



^The 5cr detection threshold was determined over a four arcsec diameter aperture using a mean sky 
noise measured in each tile. The actual depth may vary, depending on the PSF, in each tile. 
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Table 3. Summary of Optical Imaging Data 









Plate Scale 


FOV 




Field 


Telescope 


Detector 


(arcsec/pixel) 


arcmin^ 


Eandpasses 


NOAO-Cetus 


CTIO/4 m 


MosaicII 


0.27 


36x36 


V 




CFHT 


CFH12k 


0.20 


40x40 


R, L z' 


CDFS 


CTIO/4 m 


MosaicII 


0.27 


36x36 


B,V,R,I,z' 


NTT Deep . . 


CTIO/4 m 


ETC 


0.43 


29x29 


V,R,I 




CTIO/4 m 


MosaicII 


0.27 


36x36 


B,z' 


1511 


CTIO/4 m 


ETC 


0.43 


29x29 


V,R,I 




CTIO/4 m 


MosaicII 


0.27 


36x36 


B,z' 


SA22 


CTIO/4 m 


Mosaic 


0.26 


36x36 


B,V,R,z' 




CFHT 


CFH12k 


0.20 


40x40 


I 


HDFS 


CTIO/4 m 


ETC 


0.43 


29x29 


U,B,V,R,I^ 




CTIO/4 m 


MosaicII 


0.27 


36x36 


z' 



^Teplitz et al. (2001). 



This figure "f3.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0108 17 lv2 



This figure "f4.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0108 17 lv2 



Table 4: An Example of the Object Catalog of the CSDF Region 



Vk (pixels) S/G 



Tile 


ID 


X 


y 


a 


(J2000) 


(5 (J2000) 




V 


R 


/ 


z' 


H 


V 


R 


1 




z 




H 


(1) 


(2) 


(3) 


(4) 




(5) 






(6) 






(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


1 


5001 


2881 


1474 


03 


32 


30 


34 


-27 


45 


23 


70 






29 


74 










-1 


00 





00 


-1 


00 


-1 


00 


-1.00 


1 


5002 


2296 


1498 


03 


32 


30 


47 


-27 


50 


14 


99 


8 


65 


10 


71 


9 


38 









21 





52 





57 


-1 


00 


-1.00 


1 


5003 


2416 


1458 


03 


32 


30 


70 


-27 


40 


44 


92 


14 


00 


12 


99 


12 


25 


11.67 


12.91 





92 





98 





87 





45 


0.01 


1 


5004 


2151 


1510 


03 


32 


31 


15 


-27 


51 


11 


24 


8 


21 


9 


68 













27 





72 


-1 


00 


-1 


00 


-1.00 


1 


5005 


130 


1484 


03 


32 


30 


79 


-27 


46 


04 


19 


11 


63 


10 


17 


10 


45 


6.65 







12 





45 





55 





63 


-1.00 


1 


5006 


1360 


1485 


03 


32 


30 


97 


-27 


49 


18 


35 


18 


62 


13 


76 


15 


05 


15.18 


8.99 





54 





02 





02 





02 


0.01 


1 


5007 


273 


1488 


03 


32 


31 


16 


-27 


51 


23 


93 


11 


44 


11 


94 


10 


68 


10.36 


9.63 





19 





83 





79 





14 


0.24 


1 


5008 


2400 


1505 


03 


32 


30 


09 


-27 


51 


00 


03 


11 


10 


8 


22 


9 


07 









17 





13 





04 


-1 


00 


-1.00 


1 


5009 


63 


1522 


03 


32 


29 


85 


-27 


51 


06 


06 


8 


96 


10 


08 













69 





11 


-1 


00 


-1 


00 


-1.00 


1 


5010 


1209 


1506 


03 


32 


30 


03 


-27 


43 


15 


44 


12 


01 


12 


11 


11 


41 


9.24 







18 





17 





08 





53 


-1.00 


1 


5011 


484 


1514 


03 


32 


31 


21 


-27 


50 


53 


82 


10 


43 


9 


66 


10 


01 









22 





20 





19 


-1 


00 


-1.00 


1 


5012 


16 


1522 


03 


32 


30 


51 


-27 


42 


12 


76 


11 


27 


11 


94 













09 





63 


-1 


00 


-1 


00 


-1.00 


1 


5013 


105 


1469 


03 


32 


30 


07 


-27 


43 


19 


06 


30 


56 


23 


06 


16 


37 


14.21 


10.92 





04 





04 





07 





05 


0.18 


1 


5014 


82 


1457 


03 


32 


30 


96 


-27 


47 


54 


05 


23 


24 


30 


94 


16 


80 









01 





02 





10 


-1 


00 


-1.00 


1 


5015 


1838 


1471 


03 


32 


30 


59 


-27 


51 


21 


64 


17 


06 


14 


14 


11 


91 


13.99 







01 





02 





01 





01 


-1.00 


1 


5016 


128 


1525 


03 


32 


30 


63 


-27 


49 


33 


68 






7 


58 










-1 


00 





53 


-1 


00 


-1 


00 


-1.00 


1 


5017 


2072 


1495 


03 


32 


30 


98 


-27 


41 


21 


14 


10 


64 


10 


99 


16 


99 


13.42 







03 





04 





01 





53 


-1.00 


1 


5018 


1825 


1472 


03 


32 


31 


02 


-27 


42 


36 


66 


21 


70 


22 


05 


18 


76 


13.04 


12.29 





04 





04 





03 





06 


0.01 


1 


5019 


799 


1514 


03 


32 


30 


75 


-27 


43 


45 


52 


11 


39 


8 


13 


11 


09 









30 





40 





82 


-1 


00 


-1.00 


1 


5020 


24 


1493 


03 


32 


31 


28 


-27 


49 


23 


87 


16 


40 


19 


90 













04 





01 


-1 


00 


-1 


00 


-1.00 


1 


5021 


427 


1497 


03 


32 


30 


75 


-27 


43 


17 


42 


22 


85 


24 


37 


11 


38 









01 





01 





11 


-1 


00 


-1.00 


1 


5022 


2265 


1519 


03 


32 


31 


11 


-27 


47 


58 


48 






50 


51 










-1 


00 





03 


-1 


00 


-1 


00 


-1.00 


1 


5023 


1983 


1520 


03 


32 


31 


23 


-27 


45 


43 


17 


7 


00 


7 


56 


6 


85 









11 





04 





13 


-1 


00 


-1.00 


1 


5024 


1726 


1506 


03 


32 


31 


14 


-27 


42 


51 


17 


10 


64 


9 


41 


7 


46 


5.51 







91 





95 





14 





69 


-1.00 


1 


5025 


464 


1529 


03 


32 


30 


05 


-27 


39 


46 


29 


14 


98 


9 


01 













00 





10 


-1 


00 


-1 


00 


-1.00 


1 


5026 


1831 


1506 


03 


32 


31 


08 


-27 


49 


51 


10 


17 


71 


10 


29 


10 


36 


10.08 


8.47 





02 





35 





03 





02 


0.02 


1 


5027 


782 


1521 


03 


32 


30 


53 


-27 


40 


30 


33 


18 


03 


11 


39 













13 





04 


-1 


00 


-1 


00 


-1.00 


1 


5028 


1287 


1529 


03 


32 


31 


00 


-27 


40 


30 


80 


7 


80 


12 


53 


12 


34 









58 





64 





72 


-1 


00 


-1.00 


1 


5029 


1929 


1526 


03 


32 


31 


46 


-27 


50 


15 


44 


11 


60 


5 


53 













69 





68 


-1 


00 


-1 


00 


-1.00 


1 


5030 


2619 


1474 


03 


32 


30 


67 


-27 


47 


41 


70 


22 


74 


14 


05 


9 


98 









00 





44 





18 


-1 


00 


-1.00 


1 


5031 


362 


1519 


03 


32 


30 


84 


-27 


49 


31 


76 


8 


52 


8 


47 













20 





84 


-1 


00 


-1 


00 


-1.00 


1 


5032 


2455 


1497 


03 


32 


30 


21 


-27 


40 


44 


83 


16 


51 


14 


84 


15 


19 


9.25 


14.37 





01 





01 





01 





05 


0.01 


1 


5033 


2453 


1520 


03 


32 


30 


98 


-27 


43 


05 


89 


12 


53 


12 


53 


11 


71 









99 





98 





02 


-1 


00 


-1.00 


1 


5034 


271 


1537 


03 


32 


30 


92 


-27 


48 


32 


16 


11 


72 

















88 


-1 


00 


-1 


00 


-1 


00 


-1.00 


1 


5035 


845 


1500 


03 


32 


31 


37 


-27 


46 


03 


55 


22 


63 


18 


73 


9 


35 


7.56 







02 





02 





01 





15 


-1.00 


1 


5036 


434 


1507 


03 


32 


30 


98 


-27 


45 


16 


45 


11 


94 


11 


87 


17 


54 


13.30 


8.10 





56 





18 





11 





07 


0.78 


1 


5037 


2400 


1481 


03 


32 


30 


89 


-27 


42 


26 


52 


23 


00 


19 


50 


25 


50 


22.93 


13.33 





00 





00 





01 





00 


0.00 


1 


5038 


1874 


1518 


03 


32 


31 


43 


-27 


50 


42 


36 


10 


10 


8 


87 


12 


57 









36 





89 





13 


-1 


00 


-1.00 


1 


5039 


657 


1511 


03 


32 


30 


72 


-27 


46 


17 


22 


9 


38 


10 


72 


7 


46 









67 





27 





28 


-1 


00 


-1.00 


1 


5040 


1211 


1536 


03 


32 


31 


21 


-27 


43 


10 


45 






4 


58 










-1 


00 





49 


-1 


00 


-1 


00 


-1.00 


1 


5041 


1387 


1517 


03 


32 


31 


30 


-27 


48 


52 


17 


24 


95 


22 


23 


7 


49 


6.88 







03 





03 





13 





60 


-1.00 


1 


5042 


2021 


1514 


03 


32 


30 


85 


-27 


41 


57 


24 


21 


10 


10 


36 


13 


65 


19.62 







01 





04 





09 





09 


-1.00 


1 


5043 


171 


1536 


03 


32 


30 


70 


— 27 


50 


24 


73 


5 


44 

















55 


-1 


00 


-1 


00 


-1 


00 


-1.00 


1 


5044 


1160 


1503 


03 


32 


31 


21 


-27 


50 


34 


58 


14 


02 


11 


44 


12 


59 


12.11 


8.32 





02 





02 





02 





01 


0.00 


1 


5045 


1857 


1529 


03 


32 


31 


09 


-27 


45 


28 


99 


6 


34 


4 


93 


5 


25 









78 





59 





74 


-1 


00 


-1.00 


1 


5046 


582 


1531 


03 


32 


31 


18 


-27 


45 


07 


05 


9 


85 


7 


43 













55 





52 


-1 


00 


-1 


00 


-1.00 


1 


5047 


2130 


1512 


03 


32 


30 


80 


-27 


40 


04 


30 


9 


66 


10 


50 


11 


65 


7.90 


6.34 





89 





88 





70 





43 


0.33 


1 


5048 


237 


1500 


03 


32 


30 


79 


-27 


46 


37 


40 


14 


30 


12 


95 


11 


94 


9.06 







02 





03 





01 





00 


-1.00 


1 


5049 


200 


1525 


03 


32 


31 


16 


-27 


46 


56 


21 


9 


52 


8 


99 


5 


46 









30 





75 





44 


-1 


00 


-1.00 


1 


5050 


1340 


1522 


03 


32 


31 


34 


-27 


44 


46 


32 


10 


53 


10 


54 













63 





68 


-1 


00 


-1 


00 


-1.00 



49 



Table 5: An Example of the Near-Infrared and Optical Photometric Catalog of the CSDF Region'* 



Aperture (4") Best 



Tile 


ID 


X 


y 




V 








R 








1 








z' 






H 






V 








R 








1 








z' 






H 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


(20) 


(21) 


(22) 


(23) 


(24) 


1 


5001 


2881 


1474 


29 


30 






24 


65 





07 


28 


03 






23 


79 




23 


76 




29 


30 






24 


55 





07 


28 


03 






23 


79 




23 


76 




1 


5002 


2296 


1498 


24 


88 





05 


24 


22 





05 


24 


26 





09 


23 


79 




23 


76 




24 


87 





06 


24 


17 





05 


24 


16 





10 


23 


79 




23 


76 




1 


5003 


2416 


1458 


22 


25 





01 


21 


60 





00 


21 


26 





01 


20 


92 


0.02 


19 


34 


0.11 


22 


09 





01 


21 


49 





01 


21 


15 





01 


20 


82 


0.02 


19 


30 


0.11 


1 


5004 


2151 


1510 


25 


57 





11 


24 


91 





09 


27 


84 






23 


79 




23 


76 




25 


56 





10 


24 


85 





10 


27 


84 






23 


79 




23 


76 




1 


5005 


130 


1484 


24 


34 





03 


23 


51 





02 


23 


30 





04 


22 


98 


0.11 


23 


76 




24 


22 





03 


23 


41 





02 


23 


27 





05 


23 


03 


0.09 


23 


76 




1 


5006 


1360 


1485 


22 


94 





01 


23 


21 





02 


22 


33 





02 


21 


61 


0.03 


19 


52 


0.13 


22 


62 





01 


23 


03 





02 


22 


11 





02 


21 


44 


0.03 


19 


55 


0.12 


1 


5007 


273 


1488 


23 


54 





02 


22 


29 





01 


21 


54 





01 


20 


99 


0.02 


19 


53 


0.13 


23 


41 





02 


22 


18 





01 


21 


47 





01 


20 


94 


0.02 


19 


53 


0.12 


1 


5008 


2400 


1505 


24 


95 





06 


24 


60 





06 


24 


20 





09 


23 


79 




23 


76 




24 


94 





07 


24 


62 





06 


24 


24 





10 


23 


79 




23 


76 




1 


5009 


63 


1522 


25 


65 





12 


25 


07 





10 


27 


85 






23 


79 




23 


76 




25 


91 





10 


25 


30 





10 


27 


86 






23 


79 




23 


76 




1 


5010 


1209 


1506 


24 


47 





04 


23 


94 





04 


23 


73 





06 


23 


35 


0.15 


23 


76 




24 


40 





05 


23 


86 





04 


23 


64 





06 


23 


34 


0.16 


23 


76 




1 


5011 


484 


1514 


24 


58 





04 


24 


52 





06 


23 


96 





07 


23 


79 




23 


76 




24 


53 





05 


24 


53 





06 


23 


90 





08 


23 


79 




23 


76 




1 


5012 


16 


1522 


25 


10 





07 


25 


12 





11 


27 


82 






23 


79 




23 


76 




24 


90 





07 


25 


19 





09 


27 


81 






23 


79 




23 


76 




1 


5013 


105 


1469 


23 


34 





01 


22 


08 





01 


21 


02 





01 


20 


41 


0.01 


18 


85 


0.09 


22 


80 





02 


21 


69 





01 


20 


76 





01 


20 


21 


0.01 


18 


84 


0.09 


1 


5014 


82 


1457 


24 


40 





04 


23 


78 





03 


23 


68 





05 


23 


79 




23 


76 




23 


86 





04 


23 


11 





03 


23 


20 





06 


23 


79 




23 


76 




1 


5015 


1838 


1471 


24 


01 





03 


23 


36 





02 


23 


23 





04 


22 


69 


0.08 


23 


76 




23 


80 





03 


23 


13 





02 


23 


11 





04 


22 


66 


0.08 


23 


76 




1 


5016 


128 


1525 


29 


30 






25 


44 





13 


28 


03 






23 


79 




23 


76 




29 


30 






25 


29 





10 


28 


03 






23 


79 




23 


76 




1 


5017 


2072 


1495 


24 


43 





04 


24 


07 





04 


23 


62 





05 


22 


96 


0.10 


23 


76 




24 


31 





04 


24 


00 





05 


23 


03 





06 


22 


51 


0.11 


23 


76 




1 


5018 


1825 


1472 


23 


40 





02 


22 


61 





01 


21 


84 





01 


21 


31 


0.02 


19 


30 


0.11 


22 


94 





02 


22 


23 





01 


21 


52 





01 


21 


19 


0.03 


19 


24 


0.11 


1 


5019 


799 


1514 


25 


59 





10 


25 


01 





09 


24 


15 





08 


23 


79 




23 


76 




25 


32 





11 


25 


04 





09 


23 


90 





09 


23 


79 




23 


76 




1 


5020 


24 


1493 


25 


32 





09 


24 


48 





06 


27 


83 






23 


79 




23 


76 




24 


79 





09 


24 


11 





08 


27 


81 






23 


79 




23 


76 




1 


5021 


427 


1497 


24 


69 





05 


24 


38 





05 


23 


84 





06 


23 


79 




23 


76 




24 


32 





05 


23 


73 





04 


23 


93 





07 


23 


79 




23 


76 




1 


5022 


2265 


1519 


29 


30 






21 


99 





01 


28 


03 






23 


79 




23 


76 




29 


30 






20 


92 





01 


28 


03 






23 


79 




23 


76 




1 


5023 


1983 


1520 


25 


44 





09 


24 


64 





07 


24 


56 





12 


23 


79 




23 


76 




25 


35 





07 


24 


62 





06 


24 


49 





09 


23 


79 




23 


76 




1 


5024 


1726 


1506 


24 


06 





03 


23 


78 





03 


23 


73 
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'^Measurements without associated errors indicate la upper limits over a two-arcsec diameter aperture. 
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